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Optical Flow Estimation Method 

This invention relates to optical flow estimation, and, in particular, optical flow 
estimation in a compressed video data stream. 

Background of the Invention 

One of the problems of image processing lies in distinguishing foreground objects from 
background images in video data. Applications in areas as diverse as video processing, 
video compressing or machine vision rely on effective segmentation techniques to 
perform their desired tasks. Motion segmentation exploits the temporal correlation of 
consecutive video images and detects image regions with different motion. This two 
dimensional motion, usually called apparent motion or optical flow, needs to be 
recovered from image intensity and colour information in a video sequence. 

In general, depending on the target application, one can trade optimisation performance 
(accuracy) against computational load (efficiency). Some specific applications need 
very high efficiency due to real-time requirements and practical feasibility. 
Surveillance applications, such as a pedestrian detection system for underground train 
stations, are an example of a situation in which a controlled environment (fixed camera, 
controlled illumination) allied with cost requirements (large numbers of cameras and 
necessity of fast response times) is a good target for high efficiency algorithms. Such a 
system is likely to use one of the popular available video encoding standards that 
already use some form of motion estimation designed for compression purposes. 

Horn and Schunck, "Determining Optical Flow", in AL Memo 572, Massachusetts 
Institute of Technology, 1980 defines optical flow as "the distribution of apparent 
velocities of movement of brightness patterns in an image". This definition assumes 
that all changes in the image are caused by the translation of these brightness patterns, 
leading to the gradient constraint equation, involving spatial and temporal gradients and 
an optical flow velocity. 



This velocity is a two dimensional approximation of the real scene movement in the 
image plane that may be termed real velocity. The gradient constraint equation requires 
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additional constraints for resolution. Horn and Schunck (above) use a global 
smoothness term to solve this problem, while Lucas and Kanade ("An Iterative Image 
Registration Technique with an Application to Stereo Vision", Proc. Of the Imaging 
Understanding Workshop 1981 pp 121-130) use a weighted least-squares fit of local 
first-order constraints assuming that image gradient is almost constant for local 
neighbourhoods. The Lucas Kanade method generates matrix eigenvalues, the 
magnitude of the eigenvalues being directly related to the strength of edges in the image 
and the eigenvalues being used to create a confidence map of optical flow accuracy. 

A confidence map is a set of data which stores the confidence, or variance, at each pixel 
for the accuracy of the optical flow field. 

Both of the above methods are called differential methods since they use the gradient 
constraint equation directly to estimate optical flow. The largest problem of such 
differential methods is that they cannot be applied to large motions because a good 
initial value is required. 

US Patent No, 6,456,731 discloses an optical flow estimation method which 
incorporates a known hierarchically-structured Lucas Kanade method for interpolating 
optical flow between regions having different confidence values. 

MPEG-2 video encoding allows high-quality video to be encoded, transmitted and 
stored and is achieved by eliminating spatial and temporal redundancy that typically 
occurs in video streams. 

In MPEG-2 encoding, the image is divided in 16x16 areas called macroblocks, and each 
macroblock is divided into four 8x8 luminance blocks and eight, four or two 
chrominance blocks according to a selected chroma key. A discrete-cosine transform 
(DCT), an invertible discrete orthogonal transformation (see "Generic Coding of 
Moving Pictures and Associated Audio", Recommendation H.262, ISO/IEC 13818-2, 
Committee Draft MPEG-2), is applied to each 8x8 luminance block giving a matrix that 
is mostly composed of zeros (high-frequency power) and a small number of non-zero 
values. The quantization step that follows effectively controls compression ratios by 
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discarding more or less information according to the value of the quantization scale. 
Zig-zag and Huffman coding exploit the resulting high-number of zero values and 
compress the image data. 

Temporal redundancy is quite severe in video since consecutive images are very similar. 
To achieve even better compression each macroblock is compared not to its direct 
spatial equivalent in a previous image but to a translated version of it (to compensate for 
movement in the scene) that is found using a block-matching algorithm. The translation 
details are stored in a motion vector that refers to either a previous image or a following 
image depending on the picture type. 

MPEG-2 encoding defines three kinds of image data: intra-coded frame data (I pictures) 
with only spatial compression (no motion vectors), predicted frame data (P pictures) and 
bi-directionally interpolated frame data (B pictures) with motion estimation. 

I pictures only have intra-coded macroblocks (macroblocks without motion estimation) 
because they are coded without reference to other pictures. P and B pictures can also 
include inter-coded macroblocks (macroblocks where only the difference to the original 
macroblock designated by the motion vector is encoded). P pictures are coded more 
efficiently using motion compensated prediction from a past I or P picture and are 
generally used as a reference for future prediction. B pictures provide the highest 
degree of compression but require both past and future reference pictures for motion 
compensation; they are never used as references for prediction. 

Figure 1 illustrates the typical picture sequence of an MPEG-2 compressed video data 
stream. The organisation of the three picture types, I, B and P pictures, in a video 
stream is flexible, the choice being left to the encoder and being dependent on the 
requirements of the application. In the particular example shown the pictures are in the 
sequence IBBPBBP, and the arrows represent the direction in which pictures are 
estimated. I pictures are used to predict B and P pictures. P pictures are used to predict 
prior and following B pictures. B pictures are not used as references for other pictures. 
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US 6157396 discloses a system for improving the quality of digital video using a 
multitude of techniques and focussing on MPEG-2 compressed video data. The system 
aims to enhance standard compressed MPEG-2 decoding by using a number of 
additional processes, including retaining groups of pictures (GOP) and the motion 
vector information, to aid post decompression filtering in the image reconstruction (ER) 
and digital output processor (DOP). The system decompresses the image but retains the 
motion vectors for later use in the DOP. Supplemental information, such as a layered 
video stream, instructional cues and image key meta data, is used to enhance the quality 
of the decoded image through post decompression filtering. However this system relies 
on decompression of the MPEG-2 compressed video data which is disadvantageous in 
that it tends to increase computational complexity and decrease processing speed. 

It is an object of the present invention to provide fast, reasonably accurate two- 
dimensional motion estimation of a video scene for applications in which it is desired to 
avoid high computational costs and compressed digital video data is used. 

It is a further object of the present invention to provide such motion estimation which 
closely approximates the Lucas-Kanade method of optical flow estimation, but working 
only with compressed video data. 

Summary of the Invention 

According to a first aspect of the present invention there is provided an optical flow 
estimation method comprising the steps of obtaining encoded image data representative 
of an image sequence of a changing object having a motion field; extracting from said 
encoded image data first frame data blocks not incorporating motion vector encoding; 
extracting from said encoded image data second frame data blocks incorporating motion 
vector encoding; determining from said first frame data blocks confidence map data 
indicative of the edge strength within said encoded the image data and hence the 
accuracy of the motion field; deriving from said second frame data blocks smooth 
motion field data blocks in which each data block has a single motion vector and the 
magnitudes of the motion vectors are normalised; and 

updating the confidence map data on the basis of the smooth motion field data blocks to 
provide output data indicative of the optical flow of the image. 
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In one embodiment of the invention, the encoded image data is encoded in the MPEG-2 
video data format. However the invention is applicable to any compressed domain 
representation in which motion vectors are encoded. 

According to a second aspect of the present invention there is provided an optical flow 
estimation system utilising encoded image data representative of an image sequence of a 
changing object having a motion field, the system comprising first extraction means for 
extracting from said encoded image data first frame data blocks not incorporating 
motion vector encoding; second extraction means for extracting from said encoded 
image data second frame data blocks incorporating motion vector encoding; 
determination means for determining from said first frame data blocks confidence map 
data indicative of the edge strength within said encoded image data and hence the 
accuracy of the motion field; derivation means for deriving from said second frame data 
blocks smooth motion field data blocks in which each data block has a single motion 
vector and the magnitudes of the motion vectors are normalised; and updating means for 
updating said confidence map data on the basis of said smooth motion field data blocks 
to provide output data indicative of the optical flow of the image. 

According to a third aspect of the present invention there is provided computer readable 
recording medium on which is recorded an optical flow estimation program for causing 
a computer to execute the following steps: extracting, from encoded image data 
representative of an image sequence of a changing object having a motion field, first 
frame data blocks not incorporating motion vector encoding; extracting from said 
encoded image data second frame data blocks incorporating motion vector encoding; 
determining from said first frame data blocks confidence map data indicative of the 
edge strength within said encoded image data and hence the accuracy of the motion 
field; deriving from said second frame data blocks smooth motion field data blocks in 
which each data block has a single motion vector and the magnitudes of the motion 
vectors are normalised; and updating the confidence map data on the basis of the 
smooth motion field data blocks to provide output data indicative of the optical flow of 
the image. 
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Brief Description of the Drawings 

For a better understanding of the present invention and in order to show how the same 
may be carried into effect, a preferred embodiment of an optical flow estimation method 
in accordance with the present invention will now be described, by way of example, 
with reference to the accompanying drawings, in which: 

Figure 1 illustrates a typical image sequence of an MPEG-2 encoded video data stream; 

Figure 2 is a diagram providing a comparison of the Lucas Kanade method of optical 

flow estimation with that of a preferred embodiment of the present invention; 

Figure 3 illustrates the structure of the AC[1] and AC[8] coefficients within a DCT 

block of the encoded video data stream; 

Figure 4 shows a confidence map; 

Figures 5 a to 5d illustrate the steps for obtaining a smooth motion field in accordance 
with a preferred embodiment of the present invention; 

Figure 6 is a diagram illustrating the basis of a confidence map update step in such an 
embodiment; 

Figure 7 is a diagram illustrating a typical scenario where the new confidence map is a 
weighted average of the four associated DCT blocks in the I image. 
Figures 8a and 8b illustrate the effects of thresholding on a scene with no real motion; 
Figure 9 illustrates the smooth motion field of the scene of Figures 8a and 8b as 
provided by this embodiment; 

Figures 10a and 10b show the effect of noise reduction on the confidence map of the 
scene of Figures 8a and 8b; 

Figures 11a and 1 lb illustrate the effect of noise reduction on the confidence map of the 
scene of Figure 5 a; 

Figures 12a and 12b show the motion field generated by the Lucas Kanade method and 
the preferred embodiment of the present invention respectively for the scene of Figure 
5a; 

Figures 13a and 13b show the motion fields of Figures 12a and 12b without the original 
image; 

Figures 14a and 14b show further aspects of the motion fields of Figures 12a and 12b 
without the original image; 
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Figures 15a and 15b illustrate the effects of blurring caused by three confidence update 
steps applied to the motion fields of Figures 14a and 14b respectively; and 
Figure 16 is a flow diagram illustrating the steps of the method of the preferred 
embodiment of the present invention. 

Detailed Description of the Drawings 

The Lucas Kanade method of optical flow estimation noted above involves direct 
processing of pixel information. The preferred embodiment of the present invention to 
be described below closely approximates the Lucas Kanade method but working in the 
compressed video data domain, using only quantities related to the compressed video 
data. 

A parallel between the Lucas Kanade method and that of the preferred embodiment is 
drawn in Figure 2. Two tasks can be identified, namely the generation of a motion 
vector field with one vector per data macroblock and estimation of a confidence map for 
the motion vector field. Not every data macroblock has one motion vector. Some data 
macroblocks have none (intra-coded blocks) and some have two (interpolated blocks in 
B pictures). Therefore a basis is needed upon which a smooth motion field based on 
MPEG-2 motion vectors can be generated. 

The obtained motion vector field, like the optical flow field, needs some sort of 
confidence measure for each vector to be meaningful. Areas of the motion vector field 
with strong edges exhibit better correlation with real motion than textureless ones. 

A DCT block is the set of 64 DCT coefficients which result from application of a DCT 
to a data macroblock. A DCT coefficient is the amplitude of a specific cosine basis 
function, while an AC coefficient is a DCT coefficient for which the frequency in one or 
both of the dimensions is non-zero. DCT coefficients of an intra-coded data 
macroblock have a measure of edge strength in the AC coefficients. The AC[1] and 
AC[8] coefficients within such a macroblock, illustrated in Figure 3, measure the 
strength and direction of an edge within the 8x8 luminance block, and, drawing a 
parallel to the eigenvalues of the Lucas Kanade method, it is possible to use these two 
coefficients to obtain a set of confidence measures for the MPEG-2 compressed video 
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The AC[1] and AC[8] coefficients may be considered approximations to the average 
spatial gradients over a DCT data block. Let f(x,y) be the 8x8 image block 
and/ 7 (w,v) be the block's DCT. The image can be reconstructed using the inverse 
DCT: 

This reconstruction can be used for continuous values of x and y, not just integers, and 
therefore the spatial gradient f x (x,y) can be obtained by differentiating: 
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A similar expression for spatial gradient f (jc, y) can also be obtained. An average 
gradient over the complete image block can then be calculated as a weighted sum: 
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where w(x,y) is. some spatial weighting function. A popular weighting function to 
choose is a Gaussian. However, by choosing the weighting function: 

w(x 9 y) = sinl ^ — 16 - J 
for the x direction, the average gradient expression simplifies as follows: 
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where C„ and C v being functions of u and v that equal 1/V2 if u or v are 0 and C„ and 
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1 if u or v have any other value. Co and Ci are functions of u where u=0 and u=l 
respectively. 

Similar analysis can be performed for f y . Thus with appropriate weighting the AC[1] 
(F(1,0)) and AC[8] {F{§X)) coefficients can be interpreted as negative weighted 
average gradients over the DCT data block. 

To summarise the above equations, the coefficients AC[1] and AC[8] approximate to 
the average spatial gradient within a block: 

AC{\]cc-J x 

AC[S]cc-T y . 

Instead of constructing a matrix for the gradients at each pixel and then averaging, as in 
the Lucas Kanade method, the preferred embodiment of the present invention uses the 
steps of averaging the gradients (i.e. using the DCT coefficients) and then constructing a 
single matrix M' : 



This matrix will, by definition, be singular with only one significant eigenvalue: 
X x = AC[\] 2 + AC[S] 2 

which gives a measure of the confidence of the motion vector in the direction of the 
eigenvector: 



A strong eigenvalue signals a large image gradient, e.g. an edge, in this block and the 
associated eigenvector will be normal to the edge's direction. Only the magnitude of X\ 
has been used as a confidence measure, but the difference between the direction of the 
motion vector and e\ could also be analyzed. 

M r can be related to the corresponding Lucas Kanade matrix M to show that the 
resulting confidence estimates are inherently conservative. The details of this are given 
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in Appendix A. 

An MPEG-2 compressed video data stream has a variable number of DCT blocks for 
each macroblock that varies with the chroma quantization encoding parameter. Only 
the luminance DCT blocks are used in this embodiment, because there are always four 
per data macroblock (unlike the chrominance blocks that can have two, four or eight). 
Another reason is that the standard Lucas Kanade method only uses luminance values, 
so that, since the preferred embodiment attempts to approximate it, it should use only 
luminance as well. Figure 4 shows an example of a confidence map, generated from the 
AC[1] and AC[8] coefficients, of the form AC[1] 2 + AC[8] 2 , and its correlation with 
edge strength. 

As previously noted, not every data macroblock has one associated motion vector. In 
order to obtain a smooth motion field, a number of rules are implemented as follows: 

1) Macroblocks with no motion vector have the same movement as in the previous 
image. 

2) When a macroblock has two motion vectors, the one pointing back is reversed and 
added to the one pointing forward. 

3) Motion vector magnitude is normalized (motion vectors in P pictures span three 
images but this does not happen with motion vectors in B pictures, so that scaling is 
required). 

4) Skipped macroblocks in I-pictures have no movement while in P pictures they have 
movement similar to the previous block. 

Applying these rules removes the dependency on specific MPEG-2 characteristics, such 
as picture and macroblock type, and creates a motion field with one vector per 
macroblock with standardized magnitude. A spatial median filter is then applied to 
remove isolated vectors that have a low probability of reflecting real movement in the 
image. Figures 5a and 5b illustrate some stages of this process, Figure 5a showing the 
original image, and Figure 5b showing the raw vectors extracted from the data 
macroblocks of the compressed video data stream. Figure 5c shows the motion field 
after application of the above rules, and Figure 5d shows the smoothed motion field. 
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Since P and B pictures transmit mostly inter-coded blocks, the confidence map is only 
fully-defined for I pictures that typically only occur every 12 pictures. The typical 
picture sequence of Figure 1 illustrates this. A confidence map update step is required 
and the preferred embodiment exploits the fact that, even if the new exact confidence of 
a predicted data macroblock is unknown, the confidence of the original 16x16 area that 
the corresponding motion vector points at may be estimated. Figure 6 illustrates how an 
inter-coded macroblock (in this case from a P picture) relates to a corresponding I 
picture confidence map. Since motion vectors in MPEG-2 compressed video data have 
half-pixel accuracy, they rarely point to the top-left corner of a DCT data block, so that 
interpolation is required. The typical scenario will be that the new confidence is a 
weighted average of four DCT blocks. Figure 7 shows how DCT blocks may be used to 
estimate the confidence of an error block. 

There are two obvious problems that can arise from such an approach, namely excessive 
averaging and error propagation. In fact, every confidence map update step involves 
interpolating new confidences and the confidence map is only reset on I pictures that 
occur every 12 pictures in a typical IBBPBBPBBPBB sequence. In practice, this 
problem is not that serious if adequate measures are taken. An obvious measure is that 
updates should be kept to a minimum in order to avoid excessive blurring of the 
confidence map. In fact, an update is only necessary when a vector is larger than half a 
DCT data block, if a vector is smaller than this there is a low probability that that edge 
has moved significantly. Also, since motion vectors of a P picture only refer to the 
previous I or P picture, there is a maximum of three update steps, making error 
propagation less serious. Confidence map updating in B pictures may depend on the 
magnitude of motion present in the compressed video data sequence. 

Areas of the image from which edges have moved are unreliable, so that either new 
edges move to these areas or they are marked with zero confidence. This happens when 
a moving object uncovers unknown background, possibly generating random motion 
that, since it has zero confidence, is ignored. 

The method of the preferred embodiment will now be described with reference to 
Figure 16 as a sequence of steps that may be carried out within a digital processor under 
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the control of a computer program. Encoded image data representative of an image of a 
changing object, preferably in an MPEG-2 encoded format, is supplied to the digital 
processor from a video camera, for example. Step 1 corresponds to obtaining a frame 
header from the encoded image data, which, for an MPEG-2 data stream, will comprise 
information describing whether the frame is an I, B or P picture. 

At step 2 a decision is made as to whether or not the encoded video data frame 
incorporates motion vector encoding. Pictures used exclusively for reference purposes 
do not contain motion vector encoding, whereas predicted pictures contain motion 
vector information describing how to obtain the predicted picture from a reference 
picture. When the encoded video data is MPEG-2 encoded data, I pictures have no 
motion vector encoding, whereas both B and P pictures do incorporate motion vector 
encoding. The decision is therefore made on the basis of information extracted from the 
frame header. 

If the decision made at step 2 is that the encoded image data does not incorporate 
motion vector encoding then, at step 3, first frame data blocks not incorporating motion 
vector encoding are extracted. In the case of MPEG-2 encoded data, this step comprises 
extracting macroblocks of the I pictures yielding DCT coefficients. AC coefficients 
AC[1] and AC[8] are a subset of the DCT coefficients, and provide information on the 
strength and direction of edges in the real image. 

At step 4 the encoded image data not containing motion vector encoding is used to 
generate a confidence map indicative of the edge strength within the image data and 
hence the accuracy of the motion field. The AC[1] and AC[8] coefficients are used to 
create the confidence map for MPEG-2 encoded data. 

If the decision made at step 2 is that the encoded image data does incorporate motion 
vector encoding then, at step 5, second frame data blocks incorporating motion vector 
encoding are extracted. In the case of MPEG-2 encoded data, this step comprises 
extracting macroblocks of the B and P pictures. 
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At step 6, the second frame data blocks are used to derive smooth motion field data 
blocks in which each data block has a single motion vector and the magnitudes of the 
motion vectors are normalised. This is achieved by application of a set of rules which 
compensate for not every macroblock having one motion vector or motion vectors of a 
normalised magnitude. As noted above, such rules remove the dependency on the 
specific format in which the image data is encoded. 

The confidence map data determined in step 4 above is updated at step 7 on the basis of 
the smooth motion field data blocks derived in step 6. The motion vectors relate 
regions of the first frame data blocks to regions of the second frame data blocks. This 
provides a function taking the confidence map from the first frame onto the second 
frame. Because the motion vectors typically do not exactly map a DCT block in an I 
image to a DCT block in a P or B image it is necessary to interpolate across the 
confidence map as depicted in Figure 6. In one embodiment the new confidence map is 
a weighted average of the confidence values for the neighbouring 4 DCT blocks as 
shown in Figure 7. 

Steps 4 and 7 thereby together provide an updated confidence map for the smoothed 
motion vectors, leading to an estimate of the high confidence optical flow within the 
image, as performed in step 8. 

To obtain a dense optical flow field the high confidence optical flow data can be 
spatially interpolated as performed in step 9. 

Both the motion vector field and the associated confidence map are accordingly 
estimated. Magnitude and confidence map thresholding can be applied to remove the 
majority of the noisy vectors, removing glow effects of lights on shiny surfaces and 
shadows. 

Such data is somewhat sensitive to a fixed threshold and accordingly a motion 
segmentation algorithm using the motion estimation of the preferred embodiment could 
use more flexible decision methods (adaptive thresholding, multi-step thresholding, 
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probabilistic models, etc.). For the above-mentioned pedestrian detection scenario a 
fixed threshold is sufficient, and accordingly the examples referred to below use this for 
simplicity of visualization of the results. 

Using the scene of Figures 8a and 8b by way of a first example, the effect of using 
confidence thresholding to remove illumination noise and increase the accuracy of the 
motion estimation can be observed. In this specific scene, bright lights and a shiny 
pavement generate a lot of apparent motion that does not reflect real motion (Figure 9). 
The majority of the illumination noise is removed except for the noise on the handrails 
where there are strong edges and a lot of apparent motion. Figures 10a and 10b show 
the noise reduction effect of the confidence map before and after thresholding 
respectively. Such noise can be dealt with by a subsequent motion segmentation 
algorithm. A second example is shown in Figures 11a and lib with the thresholded 
low-resolution motion field (16x16 blocks) in a scene with moving objects being shown 
before and after confidence thresholding in Figures 11a and lib respectively. 

Motion fields generated according to the Lucas Kanade method provide dense motion 
maps in that a motion vector is provided for every pixel. This can be seen in Figure 12a 
where pixels with vectors pointing right are shown in white and vectors pointing left are 
show in grey for the sub-sampled Lucas Kanade motion field. Figure 12b shows the 
same scene, with the MPEG-2 smooth motion field. For ease of comparison Figure 13a 
shows this Lucas Kanade motion field subsampled to the same resolution as the 
proposed scheme, and Figure 13b shows how successfully the smooth motion field of 
the present invention approximates it. Figures 14a and 14b show both motion fields 
corresponding to another situation (that depicted in Figure 5a), these figures showing 
the subsampled Lucas Kanade and MPEG-2 smooth motion fields respectively. A close 
match may be observed between the Lucas Kanade motion field and that of the 
preferred embodiment of the present invention. 

An example of the blurring effect caused by the confidence map update step is shown in 
Figures 15a and 15b. The most severe case is presented (after three update steps) to 
show that, even with some blurring, the approximation is still quite accurate. Some 
definition is lost in the shape of the object, but most of the motion is still correct. 
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Figure 15a shows the MPEG-2 smooth motion field, and Figure 15b the Lucas Kanade 
motion field. 

The approximation to the Lucas Kanade method of optical flow estimation provided by 
the preferred embodiment of the present invention is obtained at very low 
computational cost. All of the operations are simple and are applied to a small data set 
(44*30 macroblocks, 44*30*4 DCT blocks). This, allied with minimal MPEG-2 
decoding, easily allows frame rates of over 25 images per second with unoptimised 
code. A rough comparison with Lucas Kanade method code shows that the algorithm of 
the preferred embodiment is approximately 50 times faster. With appropriate 
optimization, faster performance is possible. Given its low complexity, the method of 
the present invention may be implemented in a low cost real time optical flow estimator 
based around an inexpensive digital signal processing chipset. 

Although the smoothing effect of the weighted averaging is small, it is possible that it 
may be reduced even further. Object segmentation is likely to improve it, reducing 
blurring significantly and allowing more robust motion estimation. A more specific 
improvement would be to use low-resolution DC images for static camera applications 
where a background subtraction technique, allied with the motion estimation of the 
present invention, might be used for robust foreground object detection and tracking. 

The present invention provides a method for estimating optical flow, particularly optical 
flow as estimated by the well-known Lucas Kanade method, using only compressed 
video data. It will be appreciated by the person skilled in the art that various 
modifications may be made to the above described embodiments without departing 
from the scope of the present invention. 

Appendix A 

The Lucas Kanade and present invention matrices can be approximately characterized 
as follows: 
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and 



respectively, where E{o} represents expectation over the image pixels and the scale 
difference in M' has been ignored. It then follows from Jensen's inequality (see Cover 
and Thomas, "Elements of Information Theory" 1991, Wiley & Sons.) that for any 
direction w: 

w T M'w< w T Mw 

where (w T Mw)~ x represents the variance of the optical flow estimate in the direction w 
for the Lucas Kanade method. Thus the singular matrix M f always gives larger 
variance estimates, so that it is an inherently conservative confidence measure. 



